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Abstract. We consider a kinetic model of two species of particles interacting with a 
reservoir at fixed temperature, described by two coupled Vlasov-Fokker-Plank equations. 
We prove that in the diffusive limit the evolution is described by a macroscopic equation 
in the form of the gradient flux of the macroscopic free energy functional. Moreover, we 
study the sharp interface limit and find by formal Hilbert expansions that the interface 
motion is given in terms of a quasi stationary problem for the chemical potentials. The 
velocity of the interface is the sum of two contributions: the velocity of the Mullins-Sekerka 
motion for the difference of the chemical potentials and the velocity of a Hele-Shaw motion 
for a linear combination of the two potentials. These equations are identical to the ones 
in |OE| modelling the motion of a sharp interface for a polymer blend. 
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1. Introduction 

When a fluid mixture is suddenly quenched from a homogeneous equilibrium state into 
a thermodynamically unstable state it evolves to an equilibrium state consisting of two 
coexisting phases, each one reacher in one species. This phenomenon is called phase seg- 
regation. There are various stages during this process, starting with the formation of 
interfaces very diffuse, that sharpen with time and then move slowly driven by surface 
tension effects. In [BE LMlj the so-called late stages of the phase segregation process have 
been investigated for a kinetic model of a system of two species of particles interacting by 
a repulsive long range potential and collisions. The repulsive interaction between differ- 
ent species is modeled by a Vlasov term and the collisions by a Boltzmann Kernel. The 
evolution of the system is then ruled by two coupled Vlasov-Boltzmann equations for the 
one-particle distributions and this dynamics conserves masses momentum and energy. In 
the late stages of the coarsening process the hydrodynamics effects in this case become rel- 
evant and when the fluid is well segregated with sharp interfaces between different phases 
the interface moves in its normal direction following the incompressible velocity field solu- 
tion of the Navier-Stokes equation, while the pressure satisfy the Laplace's law relating it 
to the surface tension and curvature jBELMl] . In the present paper, we study a similar 
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kinetic model replacing the Boltzmann kernel by a Fokker-Plank operator, namely the two 
species interact with the same reservoir at fixed temperature instead of colliding each other, 
so that only the masses are conserved. We are in a situation in which the temperature 
and the momentum relax to equilibrium faster than the densities, as in the polymer blends 
where the viscosity is very large. This model can be seen as the kinetic description of a 
system of particles interacting via a weak long range potential (and with a reservoir) in 
the real space (as opposite to the lattice) or as the mean field limit of a stochastic sys- 
tem of interacting particles on the continuum (Ornstein-Ulhenbeck interacting processes). 
System of particles on the continuum are more difficult than the corresponding ones on 
the lattice because of the control of the local number of particles: the conservation law 
cannot prevent locally very high densities. Systems of this kind on the lattice have been 
introduced in a series of papers |GLPj (and references therein) to study phase separation 
in alloys and their behavior has been widely investigated. The macroscopic evolution of the 
conserved order parameter is ruled by a nonlinear nonlocal integral differential equation 
having non homogeneous stationary solutions at low temperature, corresponding to the 
presence of two different phases separated by interfaces. When the phase domains are 
very large compared to the size of the interfacial region (so-called sharp interface limit) 
the interface motion is described in terms of a Stefan-like problem or the Mullins-Sekerka 
motion depending on the time scale |GLj . In this paper we derive rigorously macroscopic 
equations for the density profiles and study by formal expansions the sharp interface limit. 
The equations for the one-particle distributions fi(x,v,r) are 

d T fi + v ■ V^/j + F{ ■ X? v fi = Lpfi i, j = 1,2, i^j, (1.1) 

where (3 is the inverse temperature of the heat reservoir modeled by the Fokker-Plank 
operator on the velocity space M. 3 

Lpfi ■= V, ■ ( M /* V *(ir))> M M = (^r 3/2 exp(-/3M 2 /2) 

and Fi are the self-consistent forces, whose potential has inverse range 7, representing the 
repulsion between particles of different species: 

F(x,r) = -V x Jdx' 1 3 U( 1 \x-x'\)jdvf J (x > ,v,T), i,j = 1,2, 1 + 3. (1.2) 

Our system is contained in a 3-dimensional torus (to avoid boundary effects) and U(r) is 
a non negative, smooth monotone function on IR + with compact support. This evolution 
conserves the total masses of the two species. Beyond the spatially constant equilibria, 
there may be other spatially non homogeneous stationary solutions. To characterize the 
stationary solutions is useful to consider the coarse-grained functional Q 

Q(Juh)'= j dxdv[(f 1 \nf 1 )(x,v) + (f 2 \nf 2 )(x,v)} + ^ J dxdv{h + f 2 )v 2 
+(3 / dxdy-y 3 U(j(x - y)) I dvfx(x,v) / dv'f 2 (y,v') 
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It is easy to see that Q is a functional decreasing in time under the Vlasov-Fokker-Plank 
dynamics. In fact, only the Fokker-Plank term gives a contribution different from zero to 
the time derivative of Q which satisfies 

|o = W wv. • (m,v„ A) m A = _ W d,d^[v„ Ap < o (1.3) 

and it is zero only and only if /j are Gaussian functions of the velocity. Hence, the 
equilibrium states are local Maxwellian with mean value u = 0, variance T = /3 -1 , and 
densities p, = J dvfi(x,v,r) satisfying 

T\og Pl {x) + J dx' 1 z U{ 1 \x-x'\)p J {x') = C i , 1 = 1,2,1^.;. (1.4) 

Moreover, the functional £ evaluated on functions of the form fi(x, v, r) = Mp(v )pi{x, r) 
with fixed total masses J dvfi, coincides apart a constant with the macroscopic free energy 
functional 

F(pi,P2) = J dx[(p l lnp 1 )(x) + (p 2 \np 2 )(x)]+P J dxdy^U(^(x - y))p 1 (x)p 2 {y) . (1.5) 

It is proved in |CCELM] . under the assumption of a monotone potential, that at low 
temperature there are non homogeneous solutions to (jl.4j) . stable in the sense that they 
minimize T . On the infinite line the non homogeneous minimizers of the excess free 
energy under fixed asymptotic values are called fronts and have monotonicity properties. 
The asymptotic values at ±00 are the values of the densities of two coexisting different 
phases at equilibrium, one reach in species 1 and the other reach in species 2. 

The macroscopic equations, which play the role of the Cahn-Hilliard equations for this 
model, are obtained in the diffusive limit: they describe the behavior of the system on 
length scales of order e^ 1 and time scales of order e~ 2 in the limit of vanishing e, where e 
is the ratio between the kinetic and the macroscopic scale. Moreover, we choose 7 = e so 
that the range of the potential is finite on the macroscopic scale. We prove in sect.s 2,3,4 
that in this limit solutions of ()1.1|) converge to solutions of the following coupled non local 
parabolic equations for the densities Pi(x,t) 

f3 2 d tPi = Ap i + f3V-{p i VU*p J ), i,j = 1,2, i^j (1.6) 

where (U * g)(x, t) = J dyU (x — y)g(y, t). These equations can be rewritten in the form of 
a gradient flux for the free energy functional T 

ft,= V.(wv£), M = r>{» »} (1.7) 

where p = {pi,p 2 ), denotes the functional derivative of T with respect to pi and M 
is the 2x2 mobility matrix. This form of the equation is very important to study the 
stability properties of the stationary solutions. Since we know that the stationary solutions 
are minimizers of the functional T , we expect to be able to prove that the system relaxes 
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to that stationary state asymptotically in time, for example using the approach developed 
in |D(JPT| for a nonconservative equation. 

To describe the late stages of the segregation process, we scale position and time as 
e _1 and e~ 3 respectively, while keeping fixed (equal to 1) 7 in the Vlasov-Fokker- Plank 
equations The width of the interface on the macroscopic scale is then of order e, 

so that in the limit e — > the interface becomes sharp. On the same scales of space and 
time the motion of the interfaces for models of alloys is given by the Mullins-Sekerka model 
|MSj . a quasi stationary boundary problem in which the mean curvature of the interface 
plays a fundamental role. We find similar results in our case, but with relevant differences. 
We choose an initial condition with an interface T and profiles for the densities given by 
front solutions with asymptotic values pf = P2 = p ± corresponding to the equilibrium 
values in the phase transition region at temperature T. In the limit e — > the difference 
of the first correction to the chemical potentials pi, ip = pi — /4 satisfies 

" A r <0(r,t) = for r G O \ T t 

. SK(r,t) 

^ t ) = p^y rGr * (1.6) 



V 



T 



2[p+-p~] 



!(p 2 -|yf)[^v r C + |^Kv r .c 



where [ ]t denotes the jump across the interface T t , p and (p are the values of total density 
and concentration respectively at equilibrium, K is the curvature in r of T t (sum of principal 
curvatures) and S the surface tension for this model (see Appendix B). This is similar to 
the Mullins-Sekerka equation but for the fact that there is an extra term determining the 
velocity proportional to v ■ V r ((r } t),r G T t where £(r, t) = (pi/4 + P2/-4~^)( r > ^ s solution 
of 

A r C(r,t) = iorreQ\T t 

[C]i = 2|^|S'A'(r,t) reT t (1.9) 

= [V ■ V r C}t 

The extra term is the normal interfacial velocity in the Hele-Shaw interface motion (jl.9|) . 
The equations for ip and ( are identical to the ones in |OE| . describing the interface motion 
of an incompressible fluid mixture driven by thermodynamic forces, modeling a polymer 
blend. A discussion on that point is in sec. 7. In sec. 6 we study the sharp interface limit 
by means of formal expansions of the Hilbert type. 

As last remark, we observe that the usual approach in literature to study the sharp 
interface limit is to start from the macroscopic equations (e.g. Cahn-Hilliard) and send to 
zero the ratio between the width of the interfacial region and the linear size of the phase 
domains. That has been studied in [MM for the macroscopic equations ()1.6JI by means of 
formal matching expansions. 



2. Macroscopic limit: expansion 



In this section we begin the study of the hydrodynamical limit for the Vlasov-Fokker- 
Plank equations (jl.lj) . We consider the diffusive scaling in which the space is scaled as e 



and time as e 2 and 7 = e so that the width of the interface is of order 1 on the macroscopic 
scale. Define 

ff(x, v, t) := /i(£ _1 x, v, e~H), i = 1,2, x G T d , v G R d . 
The equation for ff is 

dtft + \v ■ VJf + ~Ff ■ VJf = ^Lpff (2.1) 
Here Ff is the rescaled Vlasov term (jl.2)l with 7 = e: 

Ff(x,t) = -\7 X [ dx'U(\x-x'\) f dvf*(x',v,t):=-V x U*p* 

We substitute in (|2.1j) the formal power series for ff and Ff 

00 00 

/r = X> B /* (B) . ^ = E^ (n) 

n=0 n=0 

FM = -V X U* [ dvfP(x',v,t) 



We get 



^ w f) + i{ w m-„. Vl /f-^».v„/(°» 



00 ( 

E <* <wf* + " ■ ^" + " + E (**° • vjf>) - w. 



(n+2) 







n=0 k Z,Z'>0:«+i'-l=n 

At each order in e we get an equation. We write down here explicitly the first three orders: 
e~ 2 ) L p fP = 

O „ • Vjf + Ff ) ■ V w /r = hfP 

s°) ft/r + « • V,/f + if > ■ Vjf + F« ■ Vj 0) = 

From s~ 2 ) we deduce from the properties of Lp that f- ^ is the Maxwellian Mp multiplied 
by a density factor depending on x and t: 

if = Pi (x,t)Mp (2.2) 
Replacing this expression in the second equation (order e _1 ) we get 

Mpv ■ {V xPi + (3 Pi V x {U * Pj )) = Lpff ] 
So a solution has to be of the form: 

/f = Mp(A + B-v) 
where A will be fixed by the equations of the next orders and B is the vector 

B = -^{V xPi + /3piV x (U * pj)) 
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If we put these expressions of and in the e°) equation and integrate over v 
remembering that 

1 

we find the equations for the zero order densities 

1 



J 



dvMpViVj - - Uy , 



d tPl - —{A xPi + /3V x ■ { Pi V x {U-k Pj ))) = 0. (2.3) 



Now our aim is to show that a solution of equation (j2.1j) does exist and its limit as e 
goes to zero is given by (|2.2j) . with pi satisfying (|2.Hjl . We try to solve (j2.1|) in terms of a 
truncated expansion 

K 
n=0 

Replacing expression (|2.4|) for ff in equation (|2.1j) we get: 
For any n between and K 



J2 -V X U* [ dvff ] 



l,l'>0:l+l'-l=n-2 
(s) 

f i = 0, s < and for the remainder 



V v ff ] - Lpft ] = (2.5) 



dtRi + - e [v- V x Ri + Ft ■ V v Ri + B t ■ Ti] = ^L p R t - e^^A, (2.6) 



where we defined 

K 



B % = Y^e n W v ff\ T i = -V x U^ [ dvR 3 

2K-1 



n=K—l 0<1,1'<K 
l+l'—l=n 



We will find solutions to equations ()2.5|) in next section and we will study the equation 
for the remainder Ri in section 5. Here we state the results: 

Denote by (•, •)_ the following L-2 scalar product and with || • ||_ the associated norm 

(h,g)-:= / dxdvM^ 1 (v)y2[h i (x,v)g i {x,v)}. 
Jj d xm d i=12 

Put A = {Ai} i=lj2 and R = {Ri}i=x,2- Then 
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Theorem 2.1. Given a classical solution pi(x,t) of the macroscopic equations M.b]) in the 
time interval [0, T], there is a constant C depending on T, such that a unique solution to 
\2. 6)) exists and satisfies the bounds 

sup \\R(-,t)\\_ < Ce K - 1 - m \\A\\_ 

*G[0,T] 

As a consequence, 

Corollary 2.1. Under the assumptions of Theorem \2.1\ and m > 1,K — 1 — m > there 
is a positive constant e such that for e < e there is a smooth solution ff(x,v,t) to the 
rescaled Vlasov-Fokker-Plank equations \2.1\) satisfying for some constant C 

sup \\ff-Mp Pi \\_<Ce 
te[o,T] 



3. Macroscopic limit: Expansion terms 

(n) 

In this section we show existence and regularity properties of f- . For simplicity, we 
write down the proof only for K = 2, but the argument goes on for any K. The structure 
of equations (|2.5|) is very simple: they are of the form 

Lpf = h (3.1) 

with h a given function. In the Hilbert space with scalar product 

(h,g) M = dvh(-,v)g(-,v)Mp 1 
Jm d 

the kernel M = ker(Lp) is made of constants in velocity multiplied by Mp. Hence this 
equation has a solution iff h is in the orthogonal to the kernel of Lp namely iff 

J dvh(v) = 0. (3.2) 

Moreover, the solution is determined but for a term in the kernel which is of the form 
a function of x, t times the Maxwellian. Starting from the lowest order, we will see that 
h(x,v,t) = P(x,t,v)Mp with P a polynomial of the velocity with coefficients eventually 
depending on x, t. v The equation (|3.1|) can be solved uniquely in the orthogonal to the 
null space of Lp. If P is a polynomial the solution is again a polynomial of the same degree 
of P multiplied by the maxwellian Mp. In other words, if MpP G TV -1 with P a polynomial, 
then there exists a unique / 6 A/ --1 such that (J3.1|) holds. This statement can be shown by 
finding explicitly solutions to the problem 1)3.1)) for different choices of the polynomial P. 
We are interested in polynomials of degree up to the second. 

For n = 0, 1 the equations ()2.5)) are of the form ()3.1j) with h = and h = biVi respectively 
and have already been discussed in the previous section. We recall that /W can be found 
as Mp(Ai + BiVi), with Bi = — A% would be determined by the compatibility condition 
at the order n — 3. Since we are truncating the expansion at n = 2 we can safely choose 
Ai equal to zero. 
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Let us now deal with a polynomial of degree two: 

P(v) = a + biVi + CijViVj. 
By gaussian integration, the condition (|3.2|) becomes 

a + -c u = 0. (3.3) 

We look for a solution of (|3.1j) of the following type: 

f(v) = Mf3(A + B iVi + ap^). 
Plugging this ansatz in our equation we find 

d Vk (Mfsd Vk (A + + 0^^)) = M p P(v). 

Recall that d Vk Vi = 5 k i and d Vk ViVj = 5kiVj + 5kjVi] then the left hand side of the above 
equation simplifies to 

d Vk (M p (B k + C kj Vj + C ik Vi)) = 
M p (-f3v k {B k + C kj Vj + C ik Vi) + 5 kj C kj + S ik C ik ) = 
M f3 (2C ii -pB i v i -2pC ij v i v j ) 

and identifying the coefficients of the corresponding powers of v^ one gets 

Bi = Cij = ~~^jj Ci i 

with the relation a = 2Ca which is automatically verified, thanks to the compatibility 
condition ()3.3|) . In order to fix the parameter A we impose the analog of (|3.3|) : A+Ca/ (3 = 
0, namely as before we are choosing equal to zero the projection on the null space of Lp. 
Thus 

In the context of our problem the known term is always in the form of a polynomial 
multiplied by a maxwellian and the coefficients of the Vi are functions of the position. In 
the case where only first powers of v appear, i.e. P(v) = % v k , the af£ are given by 



a 



V x pi + (3piV x {U -kpj) 



here i,j = 1,2 and % ^ j. When P(v) = a® + b k vj. + c^ k VhV k the coefficients are the 
following: 

a « = dtPi + V X {U * Pj ) ■ ±{V xPi + p Pi V x (U * Pj )) 

bf = 

c hl = -^d Xh (d XkPi + (3 Pi d Xk {U * Pj )) - d Xh (U * Pj)(d Xk pi + (3 Pi d Xk {U * Pj )) 



(k) 

Summing up we denote by f i , k = 0, 1, 2 the following functions of v and x: 

= M pPi (x,t) 
if = -\M p v ■ (V xPi + (3 Pi V x (U * Pj )) 



= ~h M P d tPi + v*(u * Pj ) ■ UV xPl + f3 Pi V x (U * Pj )) 



-v ■ ±V x (v ■ (V xPi + (3 Pi (V x (U *pj))) - v- V X (U -kpj)v ■ (V x pi + (3piV x {U * pj))) 
where pi is solution of 



d tPi - ^(A xPl + pV x ■ {pN x {U * Pj ))) = 0. 



1 

The known term A\ appearing in the equation for the remainder becomes 

M = d t fP + v ■ V,/f» + ed t fP + ^ e n ~ l £ F " ] ■ V"fP 

n=l 0<l,l'<2 l+l'-l=n 

where we recall that 

F™ = -V x [ dx'U{\x-x'\) [ dv'ff 

It is easy to show that the sum over 1,1' is given by 

-V x (U* Pj )-V v fP 

indeed F^' = = F^ because the functions /f and /f belong both to A/ --1 and F^ = 
-V x {U* P j). 

In conclusion, the fj^ are always of the form Mp times a polynomial in v times a function 
of x, t which depends on the derivatives of Pi(x, t) solution of the macroscopic equations. If 
we fix an initial datum for f)1.6j) in C 2 (T d ) then the corresponding unique solution will be 
classical as shown in section 5 and the f^ n > = {f^ n '}i as well Ai will satisfy the regularity 
properties 

||/ (n) ||-<C, \\A\\_<C 



4. Macroscopic limit: Remainder 

In this section we will find a solution to equation ([2.6)1 . which is a weakly non linear 
equation ii m > 1, K — 1 — m > 0, by considering first the linear problem with the force 
term Ff assumed given so that general results will grant the existence of this linear problem 
in a suitable space. Then, a fixed point argument applies by using e as small parameter. 
From here on we will simplify notation by setting M = Mp. 

Define / = f/M and 

Lpf=±-V v -{MV V {J)). 
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Moreover, we introduce the Hilbert space associated to the L 2 scalar product (•, -) M 
weighted by the maxwellian and with || • ||m the associated norm. In this Hilbert space Lp 
is self-adjoint and non positive: 



(gi,Lf 3 g 2 ) M = (L fS g 1 ,g 2 )M 
1 

/ T d xK d M 
If we put Ri = ipiM, the equation for the remainder becomes 



{g,L p g) M 



J 

Jj d 



dxdvMg^V v ■ {MV v (g)) = -\\V v g\\ 2 M . 



M 



M 



e 2 Lpipi - e 



K-l-m 



A 

M' 



(4.1) 



In order to estimate H^Hm one multiplies the above equation by Mipi and integrates over 
x and v. So the first term on the left hand side becomes 

^II^IIm 

while the gradient with respect to the position disappear because of the periodic boundary 
conditions. 

We assume that the force terms Ff are given functions that we will call Fj and are such 
that 

I \P% I |oo — 01 pi 

Hence 



dxdv^iFi ■ V v (M^i) 



dxdvFi ■ (M^M^ipiV^i) 



< II^IIooII^IIm 



\M 



where we integrated by parts (F depends only on x) and we used Schwartz inequality. 
Now the term with the convolution of the remainder with the gradient of the potential is 
estimated in the following way: 

ixivM MA*)M {v)) ■ /™<|, - *<|) / 

dxdvpM^(v)M^(v)V v iJi(x,v) ■ J dx'dv'V x UM*{v')M*{v')i)j{x' , v') 

< \T d \ sup|p| sup|V x C/| \\ipj\\ M W^vAWm- 

As before we first integrated by parts and then we applied Schwartz inequality twice. Here 
it has been considered only the lowest order in e of the sum which constitutes Bf, the other 
terms are treated similarly. The last estimate is the one for Af 



dxdvipiAi 



d.rdrM n',— ^- < 

Ml 

2 



1 



W^Wm \\M- l A t \\ M < -{\\a\\m + Wm-'AWIj). 
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Summing up, we have 
1 
2 



K-l-m 

"^I^+IIM^II 2 ^. 



2 

Note that C\ contains powers of e greater than s~ l . Now one exploits the inequality 

£ ~ V ■(o l+ s-^<^l±^y- (4.2) 



2 v " " ~ 2 

but first we need to introduce the norm HV'Hm : = IIV'iIIm + II^Hm, so we have 

2<WH^ < -^ 2 (||V^i|lL + l|V^2||M) + (ci + £- 1 c 2 )(||V 2 || M ||V^ilUf 

+ HV'iIImIIV^IIm) + a^^dl^illMllv^ilU + II^IImIIV^Hm) 

+ —3— + + m-iMn + \\M-iA 2 \\ 2 M )} 

< {£Cl+ 2 C2) \ \\M 2 M + \\H\ 2 M ) + %\\H\ 2 m + \\H\ 2 m) 

F K-l-m 

+ -^—[(WHIm + WAW'm) + (\\M-iM\m + ||M_ a A a ||S f )] 

where we used the inequality (|4.2j) in two different ways; in fact we divided the negative 
term in two halves and then once we chose C\ = o"i and c 2 = o 2 and once we put c\ — 
and c 2 = a p. 

Multiplying by 2 both members one gets 

dt\\n 2 M<MH\\ 2 M + d 

where A = X(a p ) = a\ + (eci + c 2 ) 2 + e K ~ 1 ' m and d = e K - 1 ~ m (\\M- 1 A 1 \\ 2 M +\\M- l A 2 \\ 2 M ). 
Integrating over the time, by the Gronwall inequality: 

f{t) < Kit) + A f drf(r) < K(T) + A [ drf(r) 
Jo Jo 

fit) < K(T)e xt < K(T)e XT 

where / = HV'IIm? K(t) = Jo drd{r) is a non decreasing function of time and we used the 
initial condition /(0) = 0. 

Now consider the sequence of forces 

FW = -V S U*[ <\vYe n & ] -e m V x U* ! dvRf~ 1] 
U fo ^ 3 

with k > 1 and r\ 0) = 0. Let a k = nujx{||F :L (fc) || 0O , H^IU}, then 



a k <a + e m C J dx&v\R 



3 
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where j is chosen such that it corresponds to the maximum in the definition of a k and 



a = sup 

x€T d 



J n=0 



Write Jdv\Rf~ 1} \ = /du|M^f _1) | = f dvM*\Afty k ; using Schwartz inequality we 
get 

" dxdvlR?' 1 ^ < IT^^II^- 1 )! 



\M- 



Thus, recalling the estimate for H^IIao we can conclude that 

a k <a + £">(a fc _i) 

where the non decreasing function /j, is defined by fi(a k ) = C(\T d \K(T) exp(A(o; fc )T)) 1 / 2 . 
By induction on k we show that a k < 2a Vfc. In fact 

ol\ < a < 2a. 

Then suppose «fc_i < 2a; we have 

ttfe < a + e m /i(a fc _i) < a + e m fi(2a) < 2a 

because we applied the inductive hypothesis, exploited the monotonicity of fj, and chose e 
so small that e m ji{2a) < a. 

Denote with Sipf^ the difference — ip^ k 1 \ The equation solved by Sipf^ is 



v ■ V x (<ty< fc) ) + 



M 



+ 



M 



where is understood that Tj contains Sipf^ and no more ipj. Summing and subtracting the 



quantity ■ V v (M%p\ K 1J ) one has 
where 



(fe-i), 



(*) 



(k-i) 



-e m V x U * / du'M<ty> 



If one multiplies the equation for fo/>^ by MSip^ and integrates in space and velocities, it 
is possible to replicate the above estimates for the norm of the remainder. Only one thing 
is worth noting: the known term with Aj is now replaced by the following quantity 



dxdv84 k) SP % {k) ■ VviM^' 1 ') 



dxdvM4 k ~ 1] SP % {k) ■ v„(<^r ; ) 
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which one estimates in this way: 



dxdvM(v)(i[>f~ 1) V v 5i;f ) )(x,v) / dx'V x U{\x - x'\) j dv' M{v')5i)) ; (x>0 



<e m sup\V x U\ (^J dxdv\M^~ 1] V v c 



5^ K '\ / dx'dv' M\5ij; 



< e m \T d \ L > snp \V X U\ | |Vf _1) I \m\ | V,5# } | M _1) I 



In c the bound for ~ \\m is also present. In brief we have the following situation: 

f' k < Cf k + Qh-i 

for some C; 9 depends on e and is small as we like if m > 1, and of course f k = lllr 
with the same notation as above. By integrating in time and using Gronwall inequality we 
obtain 

h < [ T ef k ^e CT < [ T 6e CT I" Qe CT h- 2 < ... < const(6e CT T) k 



o 



thus, by a standard argument, we conclude that the sequence {ipiik)} is a Cauchy sequence 
and the limit if) is the unique solution of (j4.1|) with bounded norm HV'IIm- 



5. Limiting equation 

We follow a strategy similar to the one used in the previous section: we consider first 
a linear problem, prove existence for it and then use a fixed point argument to give the 
existence for the full non linear equation. Since we do not have at our disposal a small 
parameter we use compactness arguments and the Schauder Fixed Point Theorem [Ej . We 
seek for weak solutions in the following sense: 

Let W be the Hilbert space 

W(0, T; H\ H~ l ) :={/:/ G L 2 (0, T; H 1 ), ^ e L 2 (0, T; tf- 1 )}. 
H l (T d ) and if _1 (T d ) Sobolev spaces on the torus with norms 

I 1 2 / I 1 2 II ||2 | 1 2 i | t~7 |2 

\ V \2 = \ V \ , IMIl = Ma + I V ^2 



T d 



| \v\ |_i = sup [2(u, v) — | \u\ = dk 
ueH 1 J 

scalar product in L 2 . 

/T [H*)ll? + IKWII-i]^ 



v{k)\* 



1 + k 2 



|2 
I W 
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with v 1 = dv/dt. Let W\ be the convex subset of W 



{veW: / v(x,t) = l a.ein[0,T]} 



We say that p is a weak solution of the linear problem ()5.1|) below if for p G L 2 (T d ) and 
for all v G H 1 ^) and a.a. < t < T 

(3 2 (v, p') + (Vv, Vp + A»V[/ * A) = 

and p(-,0) = p(-). 

We remark that since p (z W implies p G C([0,T]; L 2 (T d )) we have that p(0) G L 2 (T d ). 

Theorem 5.1. For any ft G L l {T d ) and u G L 2 (T d ) there exists a unique solution in W\ 
to the following Cauchy problem 

(3 2 d t u = Au + (3V- (uV(U*h)) 

U(;0) = U(-) (5.1) 

Proof. Since h G L l (T d ) and VU as well as V 2 U are bounded we have V(Z7 * /i) and 
V 2 (U * h) in L°°([0,T] x T d ). Hence by standard arguments ;Ej there exists a solution in 
W. Since the equation is in form of divergence, the total mass is conserved so that the 
solution is in W\. 

Moreover, we have some useful a priori estimates for the solution of 15.11 (indeed the 
proof of existence can be achieved by approximation methods and these a priori estimates). 
Denote by \u\ 2 the norm in L 2 (Y d ) = J Jd dx\u\ 2 (x,t). We have that 



1 d . , 2 



2dr 



2 



1 If 

-^|Vw|l-o/ dxu(x,tWu(x,t)V(U*h)(x,t) (5.2) 



bmce h G L L (T d ) and Vf/ is bounded 



sup \V(U*h)(x,t)\ < c 

x,t 



Then, for any 5 > 



Id.., 1 o c 

< -(1 -8)^-\Vu\l + j-MuH (5.3) 



/3 2 ' " 45 
By Gronwall there exists a constant C such that 

\u\l < \u\\e ct 

so that 



[ dt\u(t)\ 2 2 < C\u\ 2 2 , [ \Vu\ 2 2 < C\u\ 2 2 
Jo Jo 
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for some constant C. Here and below C denotes a running constant. Moreover 

f 1 x — 1 

||u'||_i= sup {- / Vv[-VU * h) + —Vu} < ~| u | 2 + — \Vu\ 

veH^MvU^l JT d P P P P 



Hence 

rT 
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d*||u'(t)||ii < C\u\ 

Consider now functions u : T d — > R 2 . We define the Hilbert space W in this case as 
before, simply using as scalar product (•, •) the scalar product in L 2 (T d ; R 2 ). We use the 
same notation for W and W\. We say that p = (pi, p 2 ) is a weak solution of (jl.6|) if for all 
v e H l (T d ; R 2 ) and a.e. < t < T 

(3 2 {v, p[) + (Vv, V Pl + p Pi VU *p,) = 

and p(-,0) = p(-). 

Theorerr l5. II defines a map A from L 2 (0, T; L 2 (T d ; R 2 )) in itself by applying it to a set of 
two equations for Ui,i = 1,2 with a given term depending on gi, i = 1, 2 in the following 

way 

,3 2 c^ = Am + pVfaViUitgj)) 
Ui(-,0) = Ui(-) 
i,j = 1,2, i^j (5.4) 

We use g = {g\,g<i) and u = {u\,u<i), \g\ 2 = J2i=i 2 Ifi^li- Then, since the L 1 norm of gj 
is bounded by a constant times the L 2 norm, namely L 1 ([0, T], T d ) 6 L 2 ([0,T],T d ), there 
exists a solution u in W and we can write 

A{g) = u 

\\A(g)\\ 2 w<C\u\ 2 2 

We now prove the existence theorem for the non linear set of equations by proving that 
A is continuous and maps a closed convex set in a compact set. 

Compactness. We consider the closed and convex set X 6 L 2 ([0,T],T d ) 

X = {h : \\h\\ 2 L2{0 ^ L2) < k} 

Since A[h) is in W and W is compactly imbedded in L 2 ([0, T],T d ) the image of X is 
compact. 

Continuity. Consider g,g e L 2 (0,T; L 2 ). Let w = and u = A(g) the corresponding 
weak solutions. We have that, for i 7^ j 

(ui -Ui,u\- «•) = - — / \S7 (ui - in)\ 2 - \ ( (ui- Ui)V(ui- Ui) -V(U -k gj) 

P Jf d P JT d 

— \ UiV(ui - Ui) ■ VU * (gj - gj) (5.5) 



\~g\ u i ~ < -C\V(ui - Ui)\ 2 2 + Cijiij - Ui\\ + c 2 |^- - ffjla 
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We have used that the L 1 norm of (g — g) is bounded by the L 2 norm. Therefore, 

I \u - u\\l 2 ([o,t],l 2 ) <C\\g — g\ \l\[o,t\,l*) 

which proves the continuity of A in L 2 ([0,T], L 2 ). 

By Schauder's theorem the map A has a fixed point in L 2 ([0,T], L 2 ) which is the weak 
solution we were looking for. 

Uniqueness The proof is standard [ |GLj . 

Summarizing, we have proved the following 

Theorem 5.2. There exists a unique weak solution in W\ to the following Cauchy problem 

(3 2 d tPi = A Pi + (3V( Pi V(U* Pj )), 

ij = 1,2, i^j (5.6) 

Regularity. If VU * p G C°([0, T]; C 1 ) and p G C 2 (T d ) then the linear equation has a 
classical solution. Since the weak solution p is also in C°([0, T];L 2 ) we have that indeed 
VU -kp G C°([0, T]; C 1 ) and therefore the weak solution p corresponding to an initial datum 
in C 2 (T d ) is a classical solution. 



6. Sharp interface limit 

In this section we study the solutions of in the sharp interface limit in a 3-d torus 
Q. We introduce again the scale separation parameter e, which has the meaning of ratio 
between the kinetic and macroscopic scales. Then, we scale position and time as e~ l and 
e~ 3 , respectively, while keeping fixed (equal to 1) 7. The width of the interface on the 
macroscopic scale is then of order e, so that in the limit e — > the interface becomes 
sharp. The rescaled density distributions ff(r,v,t) = f i (e~ 1 r,v,e~ 3 t), are solutions of 

d t ft + e' 2 v ■ V r ft + e~ 2 Fl ■ V v f! = e^Lpff. (6.1) 

i^(r,t) = -V T J dr'e^Uie^lr-r'l) J dv'f*{r',v',t) =: -V r #f. 

In this section F[ depends on e through the function /? but also through the potential since 
we are keeping fixed 7. We consider a situation in which initially an interface is present. 
Since the stationary non homogeneous solutions of are given by the Maxwellian 

multiplied by the front density profiles we let our system start initially close to those 
stationary solutions and choose as initial datum ff{r,v) = Mp(v)p\, where the density 
profiles are very close to a profile such that in the bulk its values are pf, the values of the 
densities in the two pure phases at temperature T, and the interpolation between them 
on the interface is realized along the normal direction in each point by the fronts. We 
put pf = p^ and use the symmetry properties of the segregation phase transition giving 
P2 = p ± - Consider a smooth surface r C Q. Let d(r, Tq) be the signed distance of 
the point r G Q from the interface. Consider an initial profile for the densities p\ of the 
following type: at distance greater than 0[e) from the interface (in the bulk) the density 
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profiles pf(r) are almost constant equal to pf; at distance 0(e) (near the interface) we 
choose 

pt(r)=w l (e- 1 d(r,T )) + O(e) (6.2) 

where Wi(z) are the fronts, which are one dimensional solutions of (jl.4J) with asymptotic 
values pf. Since these solutions are unique up to a translation we fix a solution by imposing 
that wjo) = w 2 (0). 

Let r| be an interface at time t defined by 

r t = {reQ:pl(r,t)=pl(r,t)} 

and T be such that 1^ is regular for t G [0, T]. Let d e (r,t) be the signed distance d(r, Tf) 
of r G from the interface Tf, such that d £ > in and <i e < in f^'~, where 

= rf U U f^'~. For sake of simplicity we drop from now on the apex e. For any r 
such that |d(r, t)\ < -j^s, fc(r t ) = sup x€Ft k(x) with fc(ar) the maximum of the principal 
curvatures in x, there exists s(r) G T t such that 

i/(s(r))d(r, t) + s(r) = r 

where i/(s(r)) is the normal to the surface T t in s(r). Hence, 

is(s(r)) = V r d(r,t), reT t . 

Define the normal velocity of the interface as 

V(s(r)) = d t d(r,t). 

The curvature K (the sum of the principal curvatures) is given by K = A r d(r, t), r G Tf. 
Define, for e° small enough, 

Af(5) := {r : \d(r,t)\ < 5} 

where 5 = ±, n = max t6 [ 0iT]i0 < e < e o fc(r t ). 

We follow the approach based on the truncated Hilbert expansions introduced by Caflish 
[C] . This method, which has been used in the previous chapter to prove the hydrodynamic 
limit for the Vlasov-Fokker-Planck equation, has been improved by including boundary 
layer expansions in [ELM], to prove the hydrodynamic limit for the Boltzmann equation 
in a slab. Here we try to adapt the arguments in [ELM] to the fact that the boundary is 
not given a priori and has to be found as a result of the expansion. The Hilbert expansion 
is nothing but a power expansion in e for the solution of the kinetic equation 

oo 

/« = jV/W. (6.3) 

n=0 

Since we expect that the behavior of the solution will be different in the bulk and near the 
interface, we decompose in two parts: the bulk part p n '(r, t) and boundary terms 
which will be fast varying functions close to the interface, namely they depend on r, t in 
the following way 
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while /^(r, t) are slowly varying functions on the microscopic scale. More precisely, a fast 
varying function h(r, t) for r G Af can be represented as a function h(z, r,t), z = £ _1 <i(r, t), 
with the condition h(z,r + £z/(s(r)),t) = h(z,r,t),V£ small enough. Hence in Af we can 
write 

V r h = -vd z h + V r h\ d t h = -Vd z h + d t h; A r h = — & z h + -(V r • v)d z h + A r /i (6.4) 

£ £ £^ £ 

where the bar on the derivative operators means derivatives with respect to r, keeping 
fixed the other variables. Note that v ■ V r h(z, r, t) = 0. 

To write the expansion for the force term Ff we introduce t^ e *X]^Lo £n P^ = Y^=o £n 9i 
and Fi = — V r g- n ^. We expand also the signed distance 

oo 

d(r,t) = J^£ n d (n) (r,0 (6.5) 
i=o 

We will denote by z/ n ^ the gradient V r cr n \ with z/ := j/ -*. The condition |V r <i| 2 = 1 is 
equivalent to: 



=1 



so that g?*- ) can be interpreted as a signed distance from an interface that we denote by T t . 
As a consequence of (J6.5j) the velocity of the interface T t has the form 

oo 

We remark that giving the velocity V determines the curve evolving with it. The velocity 
V will generate an order zero interface T t . The interface generated by 'Y^ i £ l V^ 1 ' will be a 
deformation, small for small e, of T t . We define 

Af°(m) := {r : |d (0) (r,t)| < m},f t := {r : \d {0 \r,t)\ = 0},fi + H := {r : |d (0) (r, t)| > (<)0} 

and fix m so that Af°(m) C Af(S). 
We assume that in fi \ Af°(m) 



r = J2 £n f {n) - ( 6 - 6 ) 

n=0 

and that in Af°(m), the solution is of the form 

oo 

r = Y. £n f {n) ( 6 - 7 ) 



n=0 
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We will match the inner and outer expansions in z = e^ 1 m with m = £ c ,c6 (0, 1). Hence, 
we require that as z — > ±oo |CF| 

F = (F) ± + 0(e-^l) 

fV = (ffY + ^ (0) • (V^ (0) ) ± (z - df») + 0(e-^) 

f ~(2) = {f y ))± + u(0 ) . (y rf fY( z - d W) + (V r jf ^ • (-z/V) + - d«)) 
+\{drAJ? ) ) ± ^\z ~ d^)uf\z - d«) + O(e-W) 

F = + ^ " (Vr/?^^ - d«) + (V r Jf >)* • (-I/V + V V>{z - d«)) 

+^K^ fc ^ (1) ) ± ^°Vf(z - d^) 2 + (V r jf } )± • (v®(z - d«) - u^d^ - v^d®) 

+~(^aaF)>£ - rf (1) ) 3 + O(e-^l) 

where the symbol (h) for the hat functions stands for lim^ ± h(r + u£), r e f t and the 
same for the derivatives. We replace and ()6.7|) in the equations and equate terms 

of the same order in e separately in fl ± \ J\f°(m) and M°{m). We will use the notation 
Pi = I dvf!> n \ and we denote by h, h a function h(f- n ^) whenever is evaluated on f^ n \ f^. 

Outer expansion 
In Q ± \W°(m), n > 

dtft 3) +v-V r ft 1] + E fP-VJP^L^, (6-8) 

M'>0:J+/'=n-l 

with /f = 0, a < 0. 

Inner expansion 
In A/"°(m) n > we have 

M'>0:J+/'=n-2 k+k'=n 

- £ ^fV')-V,r ) + E V r ^VjT=L,/?\ (6-9) 

J,i',i">0:J+i'+i"=n M'>0:W=n-l 

with /f = 0, a < 0. 

The strategy for a rigorous proof is to construct, once the functions have been 
determined, the solution in terms of a truncated Hilbert expansion as 

N 

f = £e n / (n) +£ m R (6.10) 

n=0 

where the functions are evaluated in z = e ~ 1 d N (r,t), with d N (r,t) = £^~ 2 e n d^(r,t) and 
then write a weakly non linear equation for the remainder. In this approach it is essential 



20 



GUIDO MANZI 1 " AND ROSSANA MARRA 



to have enough smoothness for the terms of the expansion. On the contrary, they would 
be discontinuous on the border of J\f°(m) since are not exactly equal to there 
but differ for terms exponentially small in e. One can modify the expansion terms by 
interpolating in a smooth way between the outside and the inside getting smooth terms 
which do not satisfy the equations for terms exponentially small in e, that can be put 
in the remainder. With this in mind, we did not put in the equations the terms coming 
from the force such that in the convolution r is in A/"°(m) and r' in \ Af°(m). That 

is possible because the potential is of finite range. Finally, we remark that the terms f- n ^ 
of the expansion do not depend on e but for being computed on z, which depends on e 
because of the rescaling and also because the interface at time t still depends on e. The 
latter is a new feature in the framework of the Hilbert expansion due to the fact that the 
boundary is not fixed but is itself unknown. 

(n) 

In this section we show how to construct the terms fl . The argument is formal because 
we do not prove boundedness of the remainder nor the regularity properties of the terms 
of the expansion. We plan to report on that in the future. 

Now we go back to the Hilbert power series and start examining the equations order by 
order. We will find explicitly only the first three terms in the expansion to explain the 
procedure. 

Outer expansion 

At the lowest order e~ 3 (n — 0): 

L,f!° ] = 

which implies that has to be Maxwellian in velocity with variance T times a function 
pf\r,t). The latter is found by looking at the equations at the next two orders. At order 
e- 2 (n= 1): 

v-V r fV+FV.V v fV = L p fr. (6.11) 

The solution is of the form 

jf = p^Mp - M p pf\ ■ V r AS° ) (6.12) 

where /if (p £ ) = T lnpf + U 6 * p 6 j and p\ = Y^=q 5 ™/4"^ • 
The order e -1 equation (n = 2) is 

, • V r /f + ■ V,if + ■ V„jf = L p ff\ (6.13) 

The solvability condition for this equation says that the integral on the velocity of the l.h.s. 
has to be zero. By integrating over the velocity and using the explicit expression for 
we get 

-Tv r • (A (0) v r/ ,f ) = 0. 

Hence the solvability condition for the equation n = 2 gives the equation determining pf\ 
The choice of the initial data implies that the only solution of that equation is the constant 
one, with values pf in Q ± . We look at next order n = 3 to find pf^ by the solvability 
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condition. By integrating over v the equation n = 3 and taking into account that is 
Maxwellian in velocity, we get the following condition on u^ 2 \ where uf^ = J dvvf\ , 

V r -u? ] = 0. (6.14) 

Then, f- is determined, by replacing (|6.12j) in equation (|6.13j) . as 

if = -M p p®v ■ V r pf ] + pf ] M p . (6.15) 



where pf ) = Tpf ] / pf ] + <?f . 

(2) „„ K., _ + .-2 _ t-a(0)v7 ft(l) 



We use /j as given by ()6.15j) to get u 2 = —Tp\ V r p\ and plug it in ()6.14j) to get the 
equation for pf' 

A r pf ] = 0. 

We consider equation ()6.8|) for n = 3 

d t ft ) +WjP+ £ ^-VjP = L p ff\ (6.16) 

l,l>>0:l+l'=2 

whose solution is 

if = -Mf/v ■ [pf V,Af } + A (1) V,^ (1) ] + M^pfV • V r )(t; • V,)/*?* + pf ] M,. (6.17) 

The equation for pf' = Tpf* /pf' — T/2(pf^ / pf^) 2 + gP comes from the equation for 
n = 4 

l,V>0, l+l'=3 

which gives as solvability condition V r ■ if = —d t pf^ where uf^ = J dvf- 3 \ By using 
dSUH) we get 

1 V7 -C 1 ) V7 -(1) C 

Tp\ p\ p\ 

Inner expansion 
At the lowest order (n = 0) 



v.vdjr-*-Vjrd z gP = L,fP. 



In Appendix A it is proved that any solution of this equation has the form Mp(v)pf \ with 
pf^ a function of z. Plugging back in the equation we have 

d z pf ) +/3pf ) d z (U*pf)=0 8^ = 0, (6.18) 

where U is the potential U integrated over all coordinates but one. We solve this equation 
with the conditions at infinity pf, given by the matching conditions, and call Wi this front 
solution. The exponential decay of Wi has been proved for the one-component case |DOPT| 
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and the same argument should provide the proof also in this case. We can conclude that 
in Q 

ri°\r,t) = M,[w(^) Xm + (1 - Xm )pf\ 

with Xm the characteristic function of 7V°(m). This solution differs from the front solution 
Wi in Q for terms which are exponentially small in e and has the disadvantage of not being 
continuous on the border of J\f°. As explained before, it has to be modified as 

fl°\r,t) = Mp[w(^l)h(d(r,t)) + (1 - h(d(r,t)))pf\r,t)} 

with h a smooth version of Xm- 

We now find by examining the e~ 2 order (n=l) 

v ■ vdjV - v ■ Vj^dJ? - v ■ Vj^d^ = LpffK (6.19) 
The term involving i/W ■ (vd z f- ^ — V v f- ^d z g^) = (3v ■ v^Mppf^ d z pf^ = 0, because 

is solution of the lowest order equation and the bar operators vanish because pf' is 
function of z only. In Appendix A we show that the solution is necessarily Maxwellian in 



velocity so that we can write = p^Ma with p^ to be determined by the following 



equation 

d z pf ] + /3pf ] d z ~gl 1] + Pp^U * d z pf = 0. (6.20) 
Taking into account that —(3U * d z pj = d z lnwj, from the equation for the front, we get 

d z (Tpf + = dJP = 0. (6.21) 

Hence, the value of fii — p^ in z = is enough to find p^ — p^ for any z. This value 
is found as follows. From 

= T(p^)( Wi )- 1 + U*pf +K J dz'{z - z')U{z - z') Wj {z'), 

where K = A r d^ (r, t) is the curvature of the interface T t (see Appendix C) , we want to 
find ft 1 ' as determined by fip . We define the operator C as (Ch)i = Th^Wi)' 1 + U *hj. 
The previous relation reads as 

(Cp {1) )i = p\ 1] -K f dz'(z - z')U(z - z') Wj (z'). (6.22) 



The operator C has a zero mode since Cw' = 0, so that the equation (Cp^)i = hi has a 
solution only if 

/ dzhi(z)v/t(z) = 0. 

8=1,2 

The solvability condition for ()6.22|) is 

[ dzp^v/iiz) = KJ2[ dzdz'w'^iz - z')U{z - z') Wj {z'). (6.23) 

1=1,2"' 1=1,2"' 
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which implies because p\ are constant 

P?\a,r,t)[ Wl ]±™ + /4 1) (0,r,t)HlS = K(r,t) J dzdz'w>(z)(z - z')U(z - z') Wj (z'). 

In Appendix B it is shown that the sum in the right hand side is the surface tension S for 
this model, so we have (since fiui]!^ = — [1^2] 1^) 

(£« - i^mrM^itZ = K(r,t)S. (6.24) 

The matching conditions impose that p± — pip —* (M ) ± — (M f° r ^ ~~ > ±00, so that 
for rgfi 

[(M 1) ) ± -(M 1) )lM^ = ^,t)5. (6.25) 

and hence the continuity of — A2 011 the interface. 

The conservation law for the equation at the order e~ l (n = 2) will give the velocity of 
the interface. By integrating over the velocity this equation we get 

w[V + d z (v ■ up) = , (6.26) 

where the fact that and are Maxwellian in velocity is crucial for several cancella- 
tions. By integrating over z 

-V[wi]±™ = [9 ■ uftZ (6.27) 

By the matching conditions uf^ — > (u^^ at ±00, so that for r G T t 

-V[wi]±% = [v-u!?>]± ret t . (6.28) 

Summarizing what we got so far: we have constructed functions p[ harmonic in Q ± 
which satisfy (|6.25|) and (|6.28|) . For sake of convenience we will denote by p\ the func- 
tions defined in fi, not necessarily smooth, equal to (p\ ^ in Q \ Tt and such that 

lim d (o) ( J . ) t)_ + o ± fti = (Ai ) ± lr t an d the same for the derivatives. This means that Jrp 
satisfy: 

A/f> = 0, ren\t t , 

(p? - -^ ) )[p + - (T\ = K(r, t)S, r G A \ f t , 

V[ Wl ]tZ = [TpiP ■ V r ^ (1) ]i r G f t (6.29) 
where p ± = wi(±oo). Let us write the last equation as 

VP[p + - r] = l(P + ■ V r p?}± = [{(p - p)v ■ V r /4 1} ]i (6.30) 

and 

— / \ Pi(r)+p 2 (r) Pi(r)-p 2 (r) 

p( r ) = o ' ^( r ) 



2 ' T ' ' 2 
with pi{r) the step functions pi{r) := p^x + + PiX~ > the characteristic functions of 
the sets S°'(r,t) > 0, S°'(r,t) < respectively. We know, because of the symmetry of 
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the phase transition, that p is constant while (p is discontinuous in and (p(r) = ±|0| for 
r G Q ± . The previous equation implies 

2V(3[p + - p-\ = p[u ■ V r (p!p - ig>)]± + [<pv ■ VrU$> + ^)} + - (6-31) 

= p[v ■ VM ] + pV>)]± + [<pu ■ V r (p? - p?)]t (6.32) 

We introduce the function £(r, t) = (pxfli + P2/4 )( r ) t) = piflx + p 2 ) + <^(M ~ ftP) so 
that A r £(r, t) = in Q \ T t and ()6.32j) gives [u ■ V r C]l = 0. Moreover, it is discontinuous 
on Tt because of the function (p. The jump is 

C + (r,t) -C(r,t) = 2\(p\{p? -p®), r G f t 

In conclusion, ( satisfies 

A r ((r,t) = reQ\T t 
[Q± = 2\<p\SK(r,t)/[w 1 ]tZ reT t (6.33) 

= [u ■ V r C}t r G r t 

It is possible to show by using the Green identity that this problem for a given function 
K(r,t), has the unique solution 

C(r,t)= [ ds(C + -C)(s,t)v-VG(r,s) = ^^ [ dsK(s,t)u ■ VG(r, s), r G SI \ f t 

(C+ t — *) = ri4 / &>K(s, t)u ■ VG(r, 5), rGf ( 

where G is the Green function in SI. We notice that (( + + £~) = 2p(pP + pip)- 

We consider now the function £(r, £) = (pipP —p2pP){r,t) = p{pp — pP)+<p{pP +pP), 
which is discontinuous on T t and satisfies 

A r £(r,t) = refi\r 1 



[^ = y(C + + D ref, (6 . 34 ) 

T [P"Vrg]j 

2 [p+-p- 



1/ = I^A^ r e f t 



The problem is well posed because given the current configuration of the front the problem 
has a unique solution and this solution in turn determines the velocity of the front. 

In conclusion we have determined pP and pP . In Af°(m) pvp is constant equal to 
the value pP(r,t), reft, which is determined by solving the limiting equation. Hence 
pf and pi are known at this stage. As a consequence, pp are known through the 

relation pf ] = T4- + Upf in SI \ M°(m) while pP are found as solutions of (|6.22|) with 
the r.h.s. decaying to a constant as z — > ±oo and the decay is exponential if Wi do so. 
Then, a modification of the argument in |GGUlj leads to the exponential decay of p\ . 
We notice that pp is determined by (|6.22j) but for a term aw- which is in the null of 
C, with a independent of z. To fix a it is enough to put the condition pp (0, r, t)) = 
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p^\o,r,t), r G Af°(m). Since we have fixed p\ = p e 2 on T e we are allowed to choose 
p[ k) (0,r,t)) = ffi(0,r,t), r e Af°{m) for any k. 

We proceed now constructing the higher orders of the expansion. For n = 2: 

Vdjr + v ■ vd,J® + v ■ V r ff ] - d z gl% ■ V jr - d,§l% ■ V./?> 

-djV* ■ V„/?> + V r g^ • V Jf = L,f?\ (6.35) 

Again, the terms involving and are zero thanks to the previous equations. The 
matching conditions require for z large 

fi 2 \±\z\, r, t) = (ffY + v ■ (VjfYi* - d (1) ) + 0(e- Q W). 

Hence, we have to solve a stationary problem on the real line with given conditions at 
infinity. We replace in ()6.35j) = qf^ + p^Mp with J dvqf 1 = 0. This means that qf* 
is in the orthogonal to the kernel of Lp versus the scalar product 



{f,9) Mp = J dvM^fg 



We have 

Mp[Vd zWi + v ■ vd z pf + v ■ V r pf ] + Pd z gf ] v ■ v Wi + pd z gf ] v ■ vpf + pd z gfv ■ vpf 

-f3v ■ V r g < i 1) w l ] = Lpqf ] - v ■ vd z qf ] + d z gfh ■ V„# } . (6.36) 



By using the equation for the front Wi and the fact that d z p!p = together with p!p 
Tpf-lL^y + gf) weg et 



Mpv ■ v[d z pf ] + /3wid,g® + Ppf ] d z ~gf ] + ^f^] 

~(2) a ~(1) 

= Mpv ■ vwA[^- + Pg? - §(— ) 2 ] = PMpv ■ vwidj? 
Hence we can write the equation ()6.36|) in the form 

PMpu-vw&ji® = Lpqf ] -u-vd^P + dz~gt ] y-V v qf ] +(3Mpw l v-V r ~gf ) - 

- Mpv ■ V r pf ] - MpVd z Wi (6.37) 

From (j6.15)l the conditions at infinity are: 



F\±\z\r,t) = Mp [(pfT - WTv ■ (V^T + v ■ (V r pV)±(z - dP)\ + 0(e"^') 
pf } (±|^|r, t) = {pfYir, t) + v- (V r A (1) ) ± (^ - d (1) ) + 0(e- Q|21 ) 
lim / dvP(v)ql 2 \z,r,t) = -< vP(v) >p -(V r //f ^{pf ) ) ± 

z— >±oo / 
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where P(v) is a polynomial in the velocity and < - >@ are the moments of the Maxwellian 
Mp: < vP(v) >p— J dvvP(v)Mp. The matching conditions for the chemical potential 
imply that as z — > ±00 

fif = (Af°) ± + {z~ ■ (V r ^Y + O(e-^l) := (A? ) ) ± + Cf 
We have that 

^(+oo)(Af ) ) + - ^(-00)^)- = k(/2S 2) - Ci)]j£ 
where Cj = \ z< qGJ + l z> oCf . The left hand side can be written as 

dzd z \w^f ] - d)\ = J dz[ Wl d z (ti 2) - a) + <(/4 2) - a)] 

We multiply (j6.37|) by u z = v • v and integrate over v 

Widzffi = -d z [ dvv 2 z qf ] - (3v ■ uf ) (6.38) 



~(21 1 — ~f2l ~f2l 

where we used q\ G [KerLg] , z/ • V(-) = and Lpq\ = Lpf- . We have also 

Wid^fif - d) = -d z J dvv 2 z qf - (3u ■ uf - w l 8 z C l (6.39) 

By matching conditions we have that for z large 

(3v ■ uf\±\z\r,t) = -{pfYv-{VrtiY + 0{e-^) 
and that J dvv z qf vanishes at infinity. Since 

d z Q = l z<0 (u ■ Vrjlf)- + l z>0 (v ■ V r ^ (1) ) + 
we have that asymptotically in \z\ 

(59 -uf +w i d z C i = 0{e- a W) 
so that the integral over z of the r.h.s. of ()6.39|) makes sense. Then 
Wi(+oo){fif } )+ - w l {-oo)(jlf ) )- 

= -J dz[(3u ■ uf ] + wAd] + J dzw'jif - J dzw'fii (6.40) 
~ (2) 

Moreover the function v • u\ is known because 

Vw[ + d z v- uf ] = 

It follows that 



_ ~(2) Wi-p 

v ■ u\' = — 



9 p 

r,( 2 ) 



.(0) .(0) 



40) 



So only the last term involving jl\ in (|6.4(J|) needs to be computed. But we know that 

Ht ' = L pf - + # = C#> + B lfe «) + B 2 (P, (( "). (6.4i; 
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where Bi and B 2 are defined in Appendix C. The solvability condition is 

/+oo 
dz\fii ~ D iWi = 0, (6.42) 
-oo 

where Di := B\{p\ ) + B 2 (pf^). Summing on i = 1, 2 

£ 4 / = = Ei IT dzw^B^p) + fi 2 ( P 1 0) )) (6.43) 

We introduce functions pzf' defined as explained after ([6.28)1 . We have 

IplpP + P2U ( 2 ) } + - = J2 I - / ^ " ^ + + / d *<( D i ~ Ci)] = : # ( 6 - 44 ) 

We are now in position to find the first correction to the velocity of the interface, V^ 1 -*. 
This is given by the solvability condition for the boundary equation for n = 3: 

l,l'>0:l+l'=l l,l'>0:l+l'=3 
Z,Z',Z">0:Z+Z'=3 Z,Z'>0:Z+Z'=2 

After integration on t> we get 

V r • uf ] + F<9^ (1) + + 8 z {v ■ uf) + d z {y® ■ uf ] ) = 

Taking into account that for z large 

V r ■ uf ] + d z v ■ uf = (V r • ufY + v ■ (V r (P • wf ) )) ± + 0(e- a ^) 
= (V r ■ uf)* + 0(e- a|21 ) = 0(e- Q|z| ), 

(2) 

because V r ■ «, ■ = 0, we can integrate over z: 

V[P?Y™ + VU[wi]±% = ~ I dz[V r ■ uf + d z (u ■ uf)] - [i/W ■ uf (6.45) 



By the matching conditions we have as z — > ±oo 

uf~(uf)±+(z-dM)9.(V r ufY 

Let us introduce the functions 

m = l z<0 (z - d^)u • (V r (z7 • uf))~ + l z>0 (z - dF>)v ■ (V P (z? ■ uf)) + 

so that 

d z rii = l z<0 u ■ (V r (u ■ uf } ))~ + l 2>0 z> • (V r (u ■ uf)) + 

Thus 

J dz[V r ■ uf + d z {y ■ uf)] = J dz[V r ■ uf + d z ni + d z {v ■ uf - m)\ 

= J dz[V r ■ lif + d s m] + [P " uf]t =Ai + [9- uf } ]i 



A = / dz 



5 WiV r ■ V r gP + ^V r ■ V r pf } + d z rii 
J LA 

;(2) 
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— ~ (2) 

where we still have to compute Aj. Essentially we need to know V r -u\ . It can be derived 
from (|6.37|) multiplying for u± ■ v, where v± denotes one of the two directions orthogonal 
to i>, and integrating in v. After some computations we have 

Pv± ■ uf ] = -d z J dv(u± -v){v ■ v)qf ] + WiV± ■ V r g { p + -j/j. • V r pf } 

~ (2) 

In order to obtain the overlined divergence of u\ we have to sum over the two orthogonal 
directions; then remembering the asymptotic behaviour of J dvP(v)q?\ we can conclude 
that 

Here we recall that uf' depends only on quantities of the previous order in e, which are 
known. Moreover by (j6.17|) we have also 

so that 

^ • v^ 2) ]i = ntiY + v^fY - ^ ■ v r tfY- + a - i[pfV« ■ v^ (1) ]+ 

(6.46) 

We consider the functions = pijlf^ + P2pf^ an d ^ = Pin ~ P2P? ■ We have from 
ACT) 

[C (2) ]± = H(r,t) 

and by summing in equation (j6.46|) 

\v • v r c< 2 >]i = pv[p? + pZY + E • v^li + PA - [pfV 1 ) ■ v^ji} := p 

i 

Moreover, we have the identity [^Y = + ■ 

We get the velocity by taking the difference in ()6.46|) on the index i 

2[p + - P~]v m = ±[p. v r ^Y - v[Pi ] - pSY + ^[pS 1 ^ • v r ^Y- 

~[P { ^ ■ V r ^Y + A + IfpfV 1 ) ■ V r ^]±- 

-i[pf^).V r 4 1) ]l:=2[p + -p-]Q 

Notice that f r is not necessarily zero as was J Tt V. This implies that the volume 

enclosed by the interface 1$ evolving with + eV^ is not conserved. 
In conclusion, £( 2 ) and are solutions of 

A r ^(r,t) = (S 1 -S 2 )(r,t) r£0\f i 
[f (2)]+ = M(^(2)+ + C (2)-) r G f t (6.47) 

V« = Q " reft 
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and 

A r (M(r,t) = (S 1 + S 2 )(r,t) reQ\T t 

[(&}± = H(r,t) reV t (6.48) 

[u- V r C (2) ]^ = P(r,t) reT t 

Si,Q have been determined before. The terms H and P depend on $® which has been 
already found and also on S 1 ^ which is unknown and has to be determined by = dfd^ . 
These equations are different from the first order equations because the surface T t is given, 
so that we are not facing a free boundary problem. In this sense they are "linearized" even 
if the equations remain non linear. The problem is well posed because given di 1 ' on T t the 
problem has a unique solution and this solution in turn determines the velocity V^. Then, 
is found in J\f°(m) through the condition V 'd^V 'd^ = 0. The argument is analogous 
to the one in |ABCj . 

Once fip are found as solutions of these equations we can find pf* in terms of pf\ We 
have now the asymptotic values needed to solve (|6.37j) . If the solution exists, it decays 
exponentially at infinity, because the known terms have this property. This equation 
admits a solution if the conditions at infinity satisfy suitable conditions. The matching 
conditions require that the solution at infinity grows linearly. This is a problem analogous 
to the so-called Kramers problem in the half space |BCN| . It can be reduced to a Riemann 
problem with fixed asymptotic values at infinity in the following way: 
Since the solution f^' has to be approximately Af + zBf at infinity, the functions Af = 
M/ 3 ((/o? , ) ± - v ■ (pfVr/iS 1 ^) and Bf = M p {v ■ V r pf ^ have to satisfy 

v ■ v{Bf + MpftiHjfpv ■ Vrjfpf) + Mpv ■ (p^Vr^Y = Lp(Af + zBf) 
because all the other terms in (|fi.H5|) vanish at infinity. This is equivalent to 

L p {Bf) = 0, LpAf = v ■ v{Bf + Mp(3U(pf»v ■ V r pfY) + M p v ■ (pf V r //f ) ) ± 

This is true by direct inspection. Then, the problem is reduced to a well posed problem of 
finding a solution hi to eq. ()6.35|) decaying to Mp(pf )^ at infinity. 
Similar arguments lead to the computation of higher order terms. 

We conclude this section by remarking that our expansion is different from the one in 
|(J(JUlj which is more similar to a Chapmann-Enskog expansion because the terms of their 
expansion f> n ' still depend on e and are determined by equations which are nonlinear in 
the interface at every order in the sense that they are free-boundary problems determining 
for any n an interface r^ n ^ moving with velocity Yl^o^V^ ■ O ur approach is still based 
on the matching conditions and is in a way intermediate between jABC] and |Yuj . where 
it is proven the hydrodynamic limit for the Boltzmann equation in presence of shocks by 
constructing a Hilbert expansion to approximate the solution. 

7. Interface motion 

In this section we discuss the equations for the interface motion. We start by rewriting 
them in terms of the quantities ( and ip = p,^ — pip. The equation for ip is similar to 



30 



GUIDO MANZI 1 " AND ROSSANA MARRA : 



the Mullins-Sekerka equation but for the fact that there is an extra term determining the 
velocity 



A r ip(r,t) = 
,/ \ SK{r,t) 



V 



P~ 



W-\v\ 2 )[v-V r ^t + \[w-V r Q-_ 



rett\T t 

reY t 
reY t 



(7.1) 



The jump of y?z/ • V r C in the last term on the r.h.s. is indeed 2\tp\v • V r ((r, t), r E T t and 

(7.2) 



A r ((r,t)=0 

[C]± = 2\q>\SK(r,t)/[wi]l£ 



ren\r t 
rer t 
rev t 



Hence there are two contributions to the velocity of the interface: Vms, the velocity of 
an interface in the Mullins-Sekerka motion, and Vhs-, the velocity of an interface in the 
two-phases Hele-Shaw motion ()7.2|) . The latter describes the motion of a bubble of gas 
expanding into a fluid in a radial Hele-Shaw cell and is a free-boundary problem for the 
pressure P 

' A r P{r,t) = reO\T t 
[p\ + _ = CK{r,t)/{ Wl ] + _2 reT t (7.3) 



V = V ■ V r P 



Equations 17. II and 17.21 are identical to the equations in |OE] . describing the sharp interface 
motion for the dynamics of incompressible fluid mixtures driven by thermodynamic forces, 
modeling a polymer blend. In this paper the macroscopic equation is a modification of the 
Cahn-Hilliard equation for a mixture of two fluids that includes a lagrangian multiplier p 
("pressure") to take into account the constraint of constant total density 



d t Pi = V-( Pl V(i2 i +p)) 

Pl + p 2 = 1. 



1,2 



(7.4) 



This produces in the macroscopic equation for the concentration a convective term which 
in turn gives rise to the Hele-Shaw contribution Vhs to the interface motion. The macro- 
scopic equations (jl.6|) with pi = /Spi differs from the ones above for the constraint and 
hence for the pressure term. It is easy to see that the formal sharp interface limit is the 
same for both equations with V£ in the bulk a divergence-free field appearing as a velocity 
field in the equation for the total density which is constant in the bulk at the first order. 
Moreover, thermodynamic relations give that V£ = V^ 1 ) with p^ the first correction to 
the effective pressure. Hence, the role of VC is exactly the same as the lagrangian multiplier 
p in [HE] . 

We refer to |OE| for the discussion on the behavior of the interface as given in ()7.1j) 
and (|7.2|) . Here we want just to remark that the Hele-shaw motion has more conserved 
quantities than the Mullins-Sekerka motion. In fact, the former conserves the volume of 
each connected component of both phases, while the latter conserves only the total volume 
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as we can easily see by starting from 

d 



Q+\ = IV 



r 



where f2p~ is the region enclosed in the surface T. We consider now a situation in which 
there are N closed curves I\ dividing Q in N connected components f2p~ . In the Mullins- 
Sekerka problem the velocity is proportional to the jump of the normal derivative of the 
harmonic function / and this implies by using the divergence theorem 

E^rJ = £/KW]l = £/ + A /+/ A/ = 

where Q~ is the complement of Ui(f)p". U Tj). In the Hele-Shaw problem the velocity is 
proportional to the normal derivative of / and 

'''0+1= f i/-V/ = I A/ = 



dt 



In the problem 17.11 17.21 this fact has consequences on the evolution of the droplets of 
the two phases. The relative importance of the two contributions Vhs and Vms is ruled by 
the coefficients: if (p - ) -1 — << 1 (near the critical point) the Vms term dominates, 

while for deep quenches the Vhs term prevails. 



Appendix A. Uniqueness 

In this Appendix we prove that the equations for the p n \z, r, t) we examined in section 
5 have solutions whose dependence on the velocity is necessarily gaussian. We will omit 
for simplicity the dependence on the other variables. 

We consider the following set of equations for hi(z, v), i = 1, 2, (z, v ) G K x R d 

v z d z hi + Fid Vz hi = Lphi (A.l) 

where Fj = —d z f dz'U(\z — z'\)pj(z') := —d z Vi, i ^ j, Pi(z) = f dvhi(z,v), with the 
conditions at infinity 

hi(±oo, v) = M(v)pf (A.2) 

and show that it has only a solution of the form M(v)pi(z). 

Put hi = ipi(z,v)M(v)e~ l3Vi . Vi is bounded due to the assumptions on U. Then, 

Vgdgipi + Fid Vz ipi = Liipi (A.3) 

where 

Mj3 

with the conditions at infinity 

V>i(±oo,u) = e +W(±oc) p ± 
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Multiply by Mpipi and integrate over v 



^d z (v z ipi, ipi) M/3 + Fi(ipi, ^-t/)i) M(S = -{A, Li)i) M(j (A.4) 



where (h,g)M g = J dvh(v)g(v)Mp(v). We have 
Id (3 

that we write as 

d 



{v z ipi,ipi)M ~ -F^VzTp^Tp^Mp = -e ^(VuVi, ^vA)m (A.5) 



dz [Mi, ^) MfJ e- m ] = -2(V v ipi, V v iPi) M ^ Vi (A.6) 

We notice that (iv0i, -0i) (±00) = because of the boundary conditions. Hence, by inte- 
grating over z we get 

-+00 



/+00 
dze-^\\VM\l = 
-oo 



which implies V v ipi = a.e. since Vi is bounded, so that ipi = g(z), a function only of the 
position. Going back to the original equation we see that g(z) has to be the front solution. 

Next order equation. 

We discuss now equation (|6.19|) . A solution has been explicitly found as a Maxwellian 
times the density p^ l >. Suppose that there are two different solutions h and h! such that 
Ph — Ph'- Then, the equation for the difference is of the form investigated above, so that 
h — h' = 0. This means that there is a unique solution of the form Mp(v)p(z) in the class 
of solutions with fixed density p. Then, putting this expression back in the equation we 
determine p. The next order equations for have a similar form, but the solutions are 
not anymore of the form Maxwellian times a polynomial. The existence and uniqueness 
have to be proved by a different argument. 



Appendix B. Surface tension 

The surface tension for a planar interface can be defined as the difference between the 
grand canonical free energy (pressure) of an equilibrium state with the interface and a 
homogeneous one |WRj . We can call excess pressure this difference. The pressure for this 
model is 

V(n 1 ,n 2 ) = J dxp(n 1 (x),n 2 (x)) 
1 1 

p(ni,n 2 ) = T(ni logni + n 2 logn 2 ) + -^ n iU *n 2 + t^ 2 ^ * ni _ ^ lUl ~ Pi n i- 

Consider the system in a cylinder of base (2L) d ^ 1 and height M in presence of a planar 
interface dividing the cylinder in the half upper cylinder where the densities are pf, p 2 and 
the half lower cylinder with densities Pi,p 2 , where pf, pf are the equilibrium values of the 
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densities in the coexisting phases at temperature T. Then the excess pressure is given by 

m 

a = r lim f>T\d-i »i im / d ^i--- / d Vd-i dy d \p(wi,w 2 ) -p{pt,P 2 )} 

L^oo (2L) a L M^ccJ_ L J_ L J_ M 

where Wi(q) are the front solutions, smooth functions satisfying the equations 



TlogWiiq) + / dq'U{\q - g'|K(g') = Q (B.l) 
Jr 



where U(q) = L 2 dyU(^q 2 + y 2 ) and C, = /ij — T are constants determined by the condi- 
tions at infinity pf. Notice that f(pf, p 2 ) = f(PiiPt) since pf = p^ and that p\ = p 2 = p 
in the coexisting region. 

We rewrite the surface tension by using integration by part and the condition at infinity 



a = j dz[p(w u w 2 ) - p(nf,n 2 )] = - j dzz—p(w 1 (z),w 2 {z)). 



+00 



d 



00 



We have 

— p(wi,w 2 ) = T[(logWi + l)w[ + (\ogw 2 + l)w' 2 ] + -^WJ) -kw 2 

+w' 2 U -kWi + WiU * w' 2 + w 2 U -k w[] — p(w[ + w' 2 ) ■ 
By using (jB.lj) and C 1 = C 2 = C we get 

—f(wx, w 2 ) = -[—w'Jj ~kw 2 —w' 2 U ~kWi+w\U ~kw 2 +w 2 U -kw[] + (C +T)(w[+w 2 ) — p(w[+w 2 ) 
and for the surface tension, by remembering that C = p — T, 

— - J dzdz' z ^^[— w' i {z)U{z — z')wj(z') + Wi{z)U{z — z')w'j(z')] 



In conclusion, 



a = I J dzdz\z - z') Y^[w'i{z)U{z - z') Wj {z% 



Appendix C. Forces 



We show how to compute the terms gf 1 and gf 1 up to order 3. The procedure can be 
easily extended at any order. 

For a slowly varying function h(r, t) we have that 



U E *h{r,t) = I e- z U{£^\r-r'\)h(r',t)dr' 

U(\x - x'\)[h{ex',t) - h{ex,t))dx' + h(r,t) / U(\x-x'\)dx' 
U{\x-x'\) \e(x - x) • V r /i(r, t) 
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e 2 d 2 l 

+ y J2( x - x '^ x - ^jd^d~M r > *) + £4r ^ x > x ')J dx ' + 

= h(r,t)U + e 2 A r h{r,t)U + e 4 U * R h (C.l) 

where U = f U(r)dr, U — |J" r 2 U(r)dr. We have used the isotropy of [/. Hence we have 

<T = tf/f \ n = 0, 1; # = f>A (2) + ^A r pf , # = Upf + £/A r /f . 

To compute the expansion oi U £ -k h for a fast varying function r, t) it is more 
convenient to use a local system of coordinates. For a given curve T and for any point 
s 6 r we choose a reference frame centered in s with the axes 1, 2 along the directions of 
principal curvatures ki and 3 in the direction of the normal. Consider two points r and r' 
and choose the reference frame centered in s(r) : r = s(r) + ezu(r). We denote by yi and 
^ the microscopic coordinate of r and r' in this new frame. Then y\ = 3/2 = 0, — z and 
g- = £ -1 r-, the microscopic coordinates of r', are related to y[ by a linear transformation 
q[ = A ie y' e . Moreover, z' is given in terms of y[ by ( |GLj ) 

AM) =y's + E i= i )2 - 2e 2 * i y?i£] + |e 3 (E, ^V?) 2 

-^EyOfM 2 (4fci(i - <%) + 3 - 25,,)] + o(5 4 ), (c.2) 

where z' : r' = s(r') + sv(r')z' . 

We denote by hijj'^y^y'^t) the function h({eAny' e }, z'({?4}), i). We have 



(f/ £ */i)(^r,t) = / Wly-j/l)^,^,,^,*) 



/ dy'C/(|y - y'DMO, 0, ^ t) + ]- V (E/ M * 



+ iE(^^)(^r,t)) + i J2 U2 '^^V +Q 



where 



= f dyidy'MVlvs-ytf + ly'il 

Jr 2 



2 s 



We have 



/|2 ,/|2 ,/ I2)|^| 2 |^.|2_ 



<9 2 /i 2 - 2 dz' - dh dz' d 2 h dh d 2 z' 
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By using the relation between z' and y' 3 IC2I we see that the second term equals to 

dh 

— (0,0,y' 3 ,t)(ek k -e 2 2kly' 3 ). 



It is true that [B] 

fih £zi 1.2/2 K , 

l(Jl - 1 dz\z' -z)U(\z' -z\)h(z',r,t). (C.3) 



/ dy'U(\y-y'Ao,0,y' 3 ,t)J2 



To compute the contributions at different order in e we go back to the specific curve 
Tf and use the expansion d £ (r,t) = ^2 n e n d^(r,t) which implies k\ = ^2i n s n k l ( l " > and 
Afj = J2n £n ^if- m conclusion, 

(U £ *h)(z,r) = (U*h)(z,r) + e— [ dz'(z' - z)U(\z' - z\)h(z' » 

2 Jm. 

+ * D mC0 + ^ * + ^ U2,a * ^2,ii(/i) + ^C(/l)) + £ 3 5 3 

i=l,2 j^i 
3 

:= (U * r) + z nB n(h) + Rh, (C.4) 

n=l 

where i?^ is of order e 4 and 

C(h) = J dz\z' - z)U(\z> - z\)h(z',r), D 2;ij (h) = -^kf^. 
We do not write explicitly the long and uninteresting formula for B 3 . Hence we have 

~9 ( i ) ^U^pf ) +eB 1 {~pf ) ){z,r) 
§V = U*tf> + 5 1 (p1 1) )+5 2 (p1° ) ) 

~g? ) = u * pf ] + B 1 (pf ) ) + 5 2 ( P f } ) + s 3 (p1 0) ). (C.5) 
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